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1 Introduction 



Yang-Mills theory has become the candidate theory of strong interactions 
in the seventies. This confidence is based on asymptotic freedom and the 
validity of perturbation theory at short distances. The long distance/strong 
coupling regime of this theory, however, defies a sufficient understanding until 
now. More precisely, there are now several competing models well established 
on the lattice to incorporate the low energy properties of this theory. In par- 
ticular, the origin of such an distinguishing feature of strong interactions 
like quark confinement is still under debate. One of the mechanisms, the 
center vortex picture will play a prominent role in this paper. Our main 
concern, however, is the intermediate region where the operator product ex- 
pansion (OPE) is a reliable scheme to describe the onset of non-perturbative 
physics. We will point out the existence, in SU{2) lattice gluon dynamics, of 
an alternative center vortex structure which has a good chance to explain the 
nature of the gluon condensate, the most important parameter of the OPE. 

Relatively early in the history of QCD, 't Hooft has pointed out that 
choosing a particular gauge might be useful to identify the agents of confine- 
ment. In certain Abelian gauges which allow for a residual U{1) gauge 
degree of freedom, the Yang-Mills ground state appears as a dual super- 
conductor 1^ where color-magnetic monopoles of quantized charge play 
a role analogous to Cooper pairs in a superconductor. After the search for 
monopole excitations in the lattice vacua of non-Abelian gauge theory (via 
Abelian projection) has started in the late eighties 0, ^, the demonstration 
of monopole condensation, necessary to generate confinement through a dual 
Meissner effect, followed in the nineties [^, |, 0, ^ (for a caveat see [Q). This 
was demonstrated by evaluating the "disorder parameter" of confinement. 
Monopole condensation was observed for different Abelian gauges [0]. 

The idea that vortex free energies might serve as an order parameter for con- 
finement was born at the end of the seventies. It dates back to another pio- 



neering work by 't Hooft ||10[ and simultaneous work by Aharonov et al. [TT 
and was recently confirmed by a numerical investigation Yoneya |[T^ 
and Mack al. 0] were the first to construct Z{N) topological degrees of 
freedom from gauge invariant variables and pointed out that the so-called 
center vortices play an important role for the confinement of quarks. Ba- 
sically, these vortices are defined by the property that they contribute a 
non-trivial center element Za ^ 1 (among the Ai'coior^th roots of unity) to the 
Wilson loop if they are non-trivially linked to the latter. Random fiuctua- 
tions of the vortices provide the area law for the Wilson loop, the signature 
of confinement 
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A revival of the vortex picture arose with the construction of the P-vortices 
on the lattice which could be defined after choosing the so-called maximal 
center gauge (MCG) and separating out the center elements from the 
other lattice degrees of freedom, in this case by center projection. For the 
acceptance of the P-vortices as physical objects, it was essential that evi- 
dence could be presented that they are sensible degrees of freedom in the 
continuum limit |T^: the area density of the P-vortices as well as their 



binary interactions properly extrapolate to the continuum |T^. Moreover 



the P-vortex picture of the Yang-Mills ground state also gives an appeal- 
ing explanation for the finite temperature deconfinement phase transition of 



Yang-Mills theory [|T8], |l^ as the breakdown of vortex percolation. 



Subsequently, it turned out that center projection even without previous 
gauge fixing reproduces the QQ potential, also at short distances, a finding 
which has somewhat obscured the relevance of the P-vortices for forming 
the string tension [^. Moreover, in this case it was observed that other 
properties of the corresponding vortices (rather than the purely topologi- 
cal features) strongly depend on the bare lattice coupling constant (i.e. the 
lattice spacing), in a way that made it cumbersome to give a continuum 



interpretation of the apparent vortex degrees of freedom pO[. Thus, center 
projection without appropriate gauge fixing seems not to be the right way 
to avoid the apparent deficiency of the need to find an appropriate gauge, to 
fight with the corresponding technical Gribov problem etc. 

Above, vortices and the center degrees of freedom were discussed exclusively 
under the aspect of confinement. There is a general opinion that the true 
confiners are random magnetic fluxes which have some transversal extended- 
ness (thick vortices) The maximal center projection (after MCG fixing) 
however, ends up with a type of Z{2) gauge field configuration whose vortices 
(thin or P-vortices) live on the lattice scale a and are meant to localize the 
thick vortices. 

In the present paper we are going to identify another vortex structure re- 
lated to the underlying Z{2) degrees of freedom which is connected with the 
non-perturbative dynamics at short distances. In a loose sense, the 5*0(3) 
cooling method we are proposing serves to separate, among the P-vortices 
exhibited by center projection, those P-vortices which in fact represent ex- 
tended vortices of small non-Abelian action (supported by the SO (3) part of 
the gauge field) referred to as 5*0(3) vortices, from the real scale a vortices 
residing in the Z{2) part of the gauge group. Only the first ones are in the 
position to condense and are relevant for confinement. Restricting cooling to 
a finite ratio of 5*0(3) action density to the string tension a^, a residual con- 
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fining force can be defined in a lattice-scale independent way such that the 
string tension is conserved only outside a certain cooling radius and 5*0(3) 
vortices with a thickness smaller than this radius are wiped out. Infinite 
cooling leaves us with the scale a (thin) Z{2) vortex component alone. 

This seems to be a good starting point to establish a connection between 
the Z{2) gauge field content of SU{2) lattice gluon dynamics and the gluon 
condensate which (at least in the real world of SU (3) gauge theory) is an im- 
portant parameter of particle phenomenology. A large body of knowledge on 
low energy properties (resonance physics) and high energy scattering is incor- 
porated in the gluon condensate and its short range, non-local structure. In 
the case of the hadronic spectral function approaching intermediate distances 
the operator product expansion (OPE) for the current-current correla- 
tors was the first systematic framework to take into account the non-trivial 
structure of the vacuum. In this approach, non-perturbative properties of 
the Yang-Mills (or QCD) vacuum are parameterized by so-called condensates 
the values of which are fitted to the experimental hadron correlators . 

In this paper, we will demonstrate the viability of a new, gauge independent 
method to separate low from high energy 50(3) degrees of freedom and to 
suppress high energy gluons. This method exhibits a gauge independent Z{2) 
center vortex content of the SU{2) vacuum, different from the well-known 
P-vortex structure discovered by center projection. We propose that the 
gauge invariant Z(2) vortex component accounts for the average energy den- 
sity after the gluon radiation is subtracted. For this purpose, we will propose 
cooling with respect to the SO{3)=SU{2)/ Z{2) action. This kind of cooling 
eliminates the 5*0(3) part of the links which we will refer to in the following 
as gluons. We find that the SU{2) configurations, asymptotically emerging 
from the 5*0(3) cooling procedure can be considered as configurations of an 
effective Z{2) gauge model embedded into 5*f/(2). These configurations pos- 
sess an SU{2) action density which properly extrapolates to the continuum 
limit. We will argue that this action density acquires immediate importance 
as the mass dimension four condensate figuring in the OPE as gluon conden- 
sate. 

The outline of the paper is as follows. In the next section we define the 
new cooling procedure which is based on the SU{2)/Z{2) decomposition of 
the links. Section 3 contains numerical results characterizing the emerging 
gauge invariant Z{2) vortex texture. We show there that these vortices do 
not contribute to the string tension. However, they do account for the gluon 
condensate as will be shown in section 4. There we also briefly review the 
operator product approach to hadronic correlators with a special emphasis on 
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Yang-Mills condensates. The Z{2) vortex texture contribution to the gluon 
condensate is obtained there. We consider the positive plaquette model in 
section 5 focusing on the gluon condensate O4 and find that the latter is 
suppressed by one order of magnitude. Our conclusions are summarized in 
the final section 6. 



2 The gauge invariant Z(2) vortex vacuum 
texture 

2.1 Vortices and gluons 

Although the problem could be posed also for SU (3) gluon dynamics, in this 
paper, we will concentrate on 4-dimensional pure SU (2) lattice gauge theory. 
The partition function 

is a high-dimensional integral over SU (2) matrices Ui = U^,^ associated with 
the links / = {x, /i} of the lattice. The inverse bare lattice coupling 
constant [3 = 4/(yf^ [] is considered to be running with the lattice spacing 
a. Thereby, with j3 —>■ 00, convergence can be achieved towards the contin- 
uum limit for dimensionful physical quantities, such as masses, temperatures, 
string tension, condensates (all given in powers of a~^) which are extracted 
from expectation values with respect to the measure (0). We use the Wilson 
action with a density given by 

< = 1 - ^ tr f/, , (2) 

while the full action in (|l]) is a sum over all plaquettes p. Up is the usual 
ordered product Up = V Yliedp 

For the following considerations it is useful to introduce, besides the link 
variable Ux,ij, G SU{2), the adjoint link variable 

Ot = ^tr{f/.,, r'^ Ul^ r^} = O'^'IA^] , G S0{3) , (3) 

"'^The bare lattice coupling /3 should not be confused with the renormalization group 
/3-f unction (3{g). 
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which do not feel the center degrees of freedom of the hnks. These adjoint 
hnks transform under gauge transformations according 



T 



ab 



\ tr {n^r^nlr'] . 



(4) 



Comparing this transformation property with the one of continuum gluon 
fields, i.e. 



X 



(5) 



we are led to the identification of the gluon fields as the algebra fields of the 
adjoint representation, i.e. 



O 



ab 
x,ii 



exp 



{eUi{x)a} 



ab 



f^ac 



(6) 



where a is the lattice spacing. Here we here propose to distinguish between 
the gluon fields which span the SO (3) subgroup and the residual Z2 vortex 
degrees of freedom. 



2.2 Revealing the Z(2)— vortex vacuum structure 

In order to detect the inherent effective Z{2) gauge model structure, we will 
remove the gluonic (coset) degrees of freedom from SU{2) configurations 
by an appropriate cooling procedure. For this purpose, we define a gluonic 
action density per link by 




where Ox,fj.u is the plaquette calculated in terms of the 5*0(3) link elements 
Ox,tJ. (0)- The sum over z/ runs over all forward and backward directions 
orthogonal to fi. [A] is the (continuum) field strength of the (continuum) 
gluon fields A^{x) and a is the lattice spacing. 

The new cooling is performed by reducing the gluonic action, i.e. by mini- 
mization of s^l^ (1^ with respect to the fields 0^^^^. Inspired by Ref. p3|, we 
include a kind of self-restriction in our method. Further cooling of the ad- 
joint link Oxfj. is rejected iff the gluonic action is smaller than some threshold 
value 

st,<8K'a'. (8) 

Thereby k is a gauge invariant cooling scale of mass dimension one. For k = 
0, the cooling procedure completely removes the gluon fields from the SU{2) 
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lattice configurations leaving only gauge equivalents of Ox,fj. = 1- Notice that 
the standard cooling method (self-restricted or not) minimizes the full SU{2) 
action, thereby affecting center as well as coset degrees of freedom. Another 
crucial difference to the method in is that the cooling scale there measures 
the distance of the lattice configurations from classical solutions (instantons) 
while in our case the cooling scale constrains the action of the gluon fields. 

In practice, the cooling procedure works as follows in the SU{2) manifold. 
The gluonic action density s^^ (0) can be written in terms of the SU (2) 
fundamental plaquette Ux,^y, i.e. 

A local (maximal) cooling step amounts to a replacement of the link Uxfj. by 
the cooled variable U!^^^ 



(10) 

where A is a Lagrange multiplier ensuring 

UxlUl^ = 1- This local cooling 

step is disregarded iff 



1 

1 tr 

2 



< kV. (121 



The equations ( p!0| - |T^) define the cooling procedure to be applied in our in- 
vestigations reported below. Taking into account that the normalization A 
is given by A = 1/6 + O(a^), the condition ([I^ ) agrees with (|) up to order 
0{a^). One cooling sweep consists of updating once all links of the lattice 
in sequential order according (0), i.e. Ux,^ U^^^. After a finite number 



of cooling sweeps, the local constraint (|T^) is satisfied all over the lattice 
implying that there is no change in the link variables ^hj further cooling 
steps. This is how the cooling procedure stops. 



2.3 Gauge invar iance of the texture 

This cooling procedure amounts to a minimization of the S'0(3) action den- 
sity as far as tolerated by the parameter k, and brings the SU (2) plaquettes 
as close as possible to ±1. In the limit k ^ 0, this cooling eliminates the 
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^2\ p=2.3 




Figure 1: Separation of the SU{2) action into gluonic radiation (small 
and vortex vacuum texture (large /2). 



5*0(3) part of the link variables completely. Hence, the remaining field con- 
figuration can be viewed as if generated by an underlying effective Z{2) gauge 
theory. The field configurations of the latter are thin center vortices. Thus, 
the above cooling procedure extracts, in the limit k — >• 0, a structure that we 
call the Z[2) center vortex content ("vortex texture") of a given SU{2) lattice 
configuration. These vortices are given by co-closed manifolds of plaquettes 
equal to —1 and carry large Wilson action. The standard Wilson action den- 
sity can be used as a detector for the Z{2) vortex texture, even at finite k. 
Although the quantity which will be used for identifying the vortices is gauge 
invariant, one has to make sure that the same vortex structure is obtained by 
the cooling procedure of subsection when two different but gauge equiv- 
alent link configurations {f/^^}, {U^^fj.}, are analyzed. In order to see this, 
one firstly notes that the staples B^^Ofi (|TT|) transform homogeneously, i.e. 

Ux,f_i — > U^^^ : Bx^£,f^ B^,un — Bx^u^J.^x+fl ■ (13) 

Since the trace of the plaquette is gauge invariant, one finds that the cooled 
configurations obtained from {f/^^} and {Ux^fj.}, respectively, differ by the 
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same gauge transformation 

U TP Q W (14) 

The cooling procedure is thus gauge covariant and the distribution of gauge 
invariant quantities calculated on cooled configurations is independent of 
which gauge copy of the initial field configuration cooling is applied to. 

In summary, the above introduced 5'0(3) cooling facilitates a gauge invariant 
detection of the Z{2) center vortex content of (embedded in) a SU{2) lattice 
configuration. As will become clear in the following, this vortex structure 
does not coincide with the P-vortices extracted by center projecting links 
after MCG fixing. If there is coincidence with part of the latter, this part 
becomes insignificant in the continuum limit. Wc mentioned in the Intro- 
duction, that in the limit f3 ^ oo the number of (dual plaquettes forming) 
P-vortices scales like while we will find here that the corresponding num- 
ber for gauge invariant Z{2) vortices scales like a'^ (for fixed n/y/a = 0(1) 
and in the limit k — > 0). This is because our cooling method also removes 
thick vortices which have a proper support within the 5*0(3) subgroup. One 
may speculate that under cooling with respect to the 5*0(3) part of the ac- 
tion, thick vortices disappear by growing in transversal extension. Opposite 
to this, center projection (after MCG) converts thick vortices into thin ones 
which can then be detected by a large Z{2) action density at single plaque- 
ttes. This is the reason why in the following we will call thick, confining 
vortices also 50(3) vortices. 



3 Numerical results 

3.1 Clustering of SU{2) action 

Generally, the results of 50(3) cooling depend on the cooling scale k,, given 
in units of the lattice spacing. In order to express the cooling scale in physical 
units, we will relate it to the string tension, given also in lattice units. To 
get this dimensionless ratio, we will adopt the asymptotic scaling law in the 
1-loop form 

aa'{(3) = 0.12exp|-^ (/?-2.3)| , (3>2.1, (15) 

which fixes the lattice spacing a in units of the string tension a for large 
enough 
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Figure 2: Distribution of the Wilson action density on 2-dimensional hyper- 
surface after coohng steps with k = /tmax, = 620 MeV, n = 440 MeV. 



In the following, we will employ the trace of the energy momentum tensor 
6'^^ to exhibit the vortex structure of the vacuum. In the present case of a 
SU{2) Yang-Mills theory, this tensor is proportional to the Wilson action 
density (0) (see e.g. p5|)). We denote the Wilson action carried by an 
arbitrary plaquette by /2 = 1 — | tr Ux,^u- Let -P(/2) be the corresponding 
1-plaquette probability distribution. This distribution is shown in figure 
where the data come from a simulation on a 12^ lattice at (3 = 2.3. The peak 
at small values near /2 ~ can be attributed to gluon radiation, which would 
be eliminated by cooling completely only in the limit K/y/a —>■ 0. This shows 
how the constraint (p!^ ) allows to control the action contained in the form of 
gluon radiation. As expected, it decreases with stronger cooling (decreasing 
energy scale k). On the other hand, the peak at /2 ~ 2 also grows with 
increasing cooling. This shows that the contribution of the emerging gauge 
invariant Z{2) vortex texture to the total action density becomes more and 
more important with stronger and stronger cooling. As expected, the SU (2) 
action density concentrates on single plaquettes in the limit k — which are 
forming the singular vortex vacuum skeleton. Figure |^ shows the space-time 
distribution of the SU{2) action density on a 2-dimensional hypersurface of 
a generic configuration generated on a 20^ lattice at /3 = 2.3. In this plot, 
the black spots correspond to the maximum value of the action density found 
on this 2-dimensional hypersurface. One clearly observes how 5*0 (3) cooling 
leads to the clustering of action density at points where the Z{2) vortices 
pierce the considered hypersurface. 

3.2 Does the gauge invariant Z{2) vortex texture 
contribute to the string tension ? 

Let us now investigate the relation between the center vortices of the MCG 
projection (P-vortices) [|TB| and the gauge invariant Z{2) vortex texture 
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Figure 3: The static quark QQ force as function of the distance r between 
quark and anti-quark for full non-Abelian and 5*0 (3) cooled configurations, 
at different k/a/o" checked with respect to scaling for various (3. 



which is defined by the Z{2) gauge fields remaining after the 5*0 (3) cool- 
ing described above. The contribution of these Z{2) vortices to the string 
tension is useful to look at. It was established that configurations the links of 
which were projected onto center elements after MCG fixing - the P-vortex 
configurations - reproduce the full string tension a to good accuracy (for 
most recent results, which also cover a discussion of the practical Gribov 



problem see [0, |2 



For comparison, we investigate the static QQ attractive force and contrast 
the result from cooled configurations (at different cooling scales representing 
the remaining 50(3) action) to the result obtained in full SU{2) Yang-Mills 
theory. In each case, we used /5- values ranging from 2.1 to 2.6 to check for 
scaling of the force as a function of the distance. Our results for a 12^ lattice 
are shown in figure We find that the data points (for no cooling and a 
given cooling scale) obtained at different /3 fall on top of the same curve, re- 
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spectively, thus establishing proper scahng also for the potential gained from 
cooled configurations. This shows the advantage to have a renormalization 
group invariant formulation for restricted cooling. We find that the cooling 
procedure strongly affects the force in a range of distances growing with de- 
creasing cooling scale nj^. At short distance this was expected since the 
behavior at small r is dominated by the exchange of gluons, which is already 
partially (cf. figure [l|) eliminated by cooling. Moreover, for stronger cooling 
(smaller cooling scale k) the value of the full string tension is approached only 
asymptotically, such that lowering the cooling scale k/ y/a shifts the asymp- 
totic region to larger distances r. The explanation is that S'0(3) cooling 
washes out the 50(3) vortices. We conclude that the gauge invariant Z(2) 
vortex texture (which survives in the limit k ^ of unrestricted cooling) is 
not related to the confinement property of the SU{2) vacuum. In particu- 
lar, it is not identical with the P-vortex ensembles constructed by means of 
center projection. 

After the major part of the confining vortices is removed by cooling, it re- 
mains to be clarified which specific physical significance could be assigned 
to the exposed Z{2) vortex ensembles. Following the original version of this 



paper, it was observed in that the masses of the low-lying glueballs, 
and 2+, are rather insensitive against the here proposed cooling. Below, we 
will give further arguments supporting the idea that these configurations are 
not just lattice artefacts but play an important role in the non-perturbative 
physics at intermediate distances dealt with in the OPE. 



4 Hadronic correlation functions and OPE 
4.1 The operator product expansion 

In this section, we will briefly review the operator product expansion (OPE), 
which is the framework to discuss the properties of hadronic resonances p2| , 
and the inherent ambiguities. To be specific, consider the hadronic current 
correlation function in the vector meson channel, 

M'^'iq) = I j d^x e*^^ {n\ T f{x) j'^(O) , (16) 

where \VL) denotes the true ground state of Yang-Mills theory. Current con- 
servation implies 

M^^(g) = (gV - q^g^^) M{q') . (17) 
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vector meson 




Figure 4: Two Feynman diagrams contributing to the current-current corre- 
lation function in the vector meson channel; radiative gluon (a) and conden- 
sate contribution (b). 



According to Wilson [28|, the bilocal, time-ordered product of two operators 



at short distance can be written in terms of local operators 

oo 

Tf{x)f{0) = Y^CM 6„(0) , do = 1 



n=0 



where n labels the canonical mass dimension of the local operators 0„. The 
Wilson coefficients C„(x) contain the singular behavior occurring when the 
point splitting is removed as x — > 0. The vacuum expectation values of the 
operators 0„ correspond to physical observables, which are called conden- 
sates. Hence, sandwiching (plSf ) with an arbitrary trial state {ip) and resorting 
to a dimensional analysis, one finds 

M{q') = Ctiq') + ^ + ^ + {l/q') , Ot := (^|On|V^) • 

^ ^ (19) 

Since an infinite number of degrees of freedom is encoded in the wave function 
\ip), the matrix element of the unit operator, = is generically 

depending on the state under discussion. Usually this is included in the 
coefficient function. Therefore the function CQ^q^) in (|T^ does depend on 



the state (see below and p9l for an illustration). This fact is just a refiection 
of the scale anomaly I^Oj. The other Wilson coefficient functions are defined 
in ([T8|) without reference to a particular state. Note also that a logarithmic 
dependence of Co'(g^) on is not excluded by the dimensional analysis. In 
particular, choosing \ip) = |0), the perturbative vacuum satisfying 0^ = 
for n > 0, yields 

C'.iq') = Mpe.(g2) , (20) 
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where Mpert is the correlator ([T7| ) calculated by summing perturbative dia- 
grams up to a finite order using tree-level propagators and vertices. In this 
case, one finds up to second order in the gauge coupling constant g 



C„V) = -5^5(1 + ^^)14 + ("3, (21) 
where u is the renormalization point and as = (7^/47r. 

It has been known for some time that the condensates are ambiguous since 
a splitting like (|19D into a so-called perturbative part and contributions 
from the OPE corrections is not well-defined (see e.g. [PT|). This conclu- 



sion is based on the observation that a partial re-summation of perturba- 
tive diagrams to all orders can generate, after renormalization, contributions 
exp(— p-) to correlators. Containing essential singularities in the coupling 
constant, these contributions are imitating OPE condensate corrections (see 
e.g. [0). These results do not invalidate the considerations above since the 
re-summation of perturbative diagrams to all orders bears the potential to 
change the properties of the reference state under discussion. This has 



been illustrated for the two-dimensional Gross-Neveu model in Ref . p9[ . The 
model is solvable in the large Nf limit where Nf is the number of fermion 
flavors. In leading order of the l/A^/-expansion, the ladder re-summation 
of perturbative diagrams results in an exact, non-perturbative gap equa- 
tion which admits a non-trivial solution corresponding to the true ground 
state for Nf 00, which we call \BCS). Doing this re-summation of all 
perturbative diagrams which contribute in the large Nf limit yields 

M{q') = C^^^q^) + ^ + ^ + O (1/g^) . (22) 

In this example, the numbers O^'"^ , which vanish in the infinite Nf limit, 
account for correlations which are not present in the wave function \BCS) 
and which become non-negligible at finite Nf. Moreover, C(f'"'^(g^) con- 
tains already the power-like OPE corrections. If one does not perform the 
summation over an infinite number of diagrams, one obtains from the OPE 

Miq^) = C',{q^) + % + ^ + (1/g^) , (23) 

where Cg(g^) is calculated from a large but finite number of perturbative 
diagrams (see (PO])). The superscript signals that the condensates 0° carry 
all the information on the true vacuum state, since contributions which could 
imitate a condensate are lacking in CQ(g^). In particular, the parameter O2 
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reflects the dynamical generated fermion mass which is present in the true 
vacuum of the infinite Nj limit but goes beyond perturbation theory. 

Both (p3D and (p2D, respectively, represent the unique answer for M(g2) at 
large values of g^. The example of the Gross-Neveu model demonstrates 
that choosing for {ip) the true ground state \Q) does not necessarily imply a 
unique definition of the condensates O^. In addition, usually one demands 
that, with Mper given in 

C^{q') = Mper(g^) . (24) 

which is indeed satisfied in the case of (^). The summation of a finite (even 
large) number of perturbative diagrams cannot imitate the renormalization 
group invariant dependence of the condensates on the coupling g. Therefore, 
the choice (0) is specific in as far it shifts the maximum content of informa- 
tion on the non-trivial properties of the true vacuum \Q) to the condensates 



4.2 The Yang— Mills condensate on the lattice 

The local condensates On entering the OPE ([18|) are required to be physical 
observables. Hence, in the case of Yang-Mills theory, they consist of gauge 
invariant and renormalization group invariant quantities. In the present pa- 
per, we will concentrate onto the condensate O4 of canonical mass dimension 
four. 

Using the trace of the energy momentum tensor 6^^ for the operator O4 meets 
the gauge invariance condition. However, it has been known for a long time 
that the contribution from gluon radiation yields divergent results for the ex- 
pectation value (O'^ij) , hence, violating the second criterion of renormalization 
group invariance. For later convenience, we illustrate this observation for the 
case of lattice regularization. Calculating 6*^^ by studying scale variations, 
one finds 

("",.) = A (1 - 5 f,) . (25) 

^ a* a ma 2 

where /5(a) is the running bare lattice coupling. Using (lattice) perturbation 
theory one finds that 

(1 - ^ tr = A + 0(l//?2) . 

Given (0), this implies that (0^^) strongly diverges in the continuum limit 
a ^ 0. 
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In order to arrive at a sensible (i.e. gauge and RG invariant) definition of 
O4, one closely follows the perturbative approach augmented by OPE cor- 
rections and defines 



:= (: e% ■■) ■■= (^^)sub (26) 

where (^'^^)sub denotes the expectation value from which the contribution 
from gluon radiation is subtracted. A technique subtracting perturbative 
contributions was firstly employed in lattice calculations in the early eight- 
p2| , p3| , 0, and was pursued to high order in a recent study of the OPE 



les 



corrections 



By subtracting perturbative gluonic parts, both, the continuum approach 
which distinguishes between perturbative gluon radiation and non-perturbative 
vacuum properties as well as the lattice approach, have the shortcoming that 
they do not specify the remnants of the vacuum which form the gluon con- 
densate once gluonic contributions have been subtracted to all orders. 

Here, we define the condensate through the subtraction of the contribution 
of the gluons (i.e. subtracting the 5*0(3) coset fields) by means of the above 
cooling procedure rather than by subtracting perturbation theory, i.e., 

O4 := lim04(K), (27) 

24 / 1 \ 
0,{K)a^ = (1 _ trf/p) 

Vr"^ \ 2 / 50(3) cooled with scale k 

This definition amounts to the choice - in the sense of eq. (|^) - 

C^{q') = Mso(3)(g'). (28) 

We point out that Mso(3) not necessarily contains only perturbative gluon 
contributions. Hence, in general, Mso(3) 7^ Mp^r- 

This identification offers the appealing feature that the condensate of mass 
dimension four, O4, could be given by the vacuum energy density of the Z2 
gauge system which remains after suppressing the contribution of the coset 
fields by the 5*0(3) cooling method. It remains to show that the condensate 
O4 (^) properly extrapolates to the continuum limit a — 0, which will be 
done in the next subsection. 



4.3 Numerical results for the Z{2) vortex induced 
gluon condensate 

In section 0, we have demonstrated a string like clustering of SU (2) action 
density which corresponds to a gauge invariant Z{2) vortex vacuum texture. 
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Figure 5: O4 in units of the lattice spacing a extracted from 50(3) cooled 
configuration as function of f3. 



In this subsection, we will show that, in addition to gauge invariance, this 
action density, which is carried exclusively by the Z{2) gauge fields, possesses 
the correct renormalization group dependence on the SU{2) inverse coupling 
(3 inherited from the full SU{2) lattice gluon dynamics. Hence, this action 
density can be considered as a physical observable and can be used to define 
the mass dimension four condensate O4 as outlined in the previous subsection. 

Recalling (^Sf ), O4 should be compared with what is defined in the literature 
as the gluon condensate 



4^2 ■ tJ-i^ 



(29) 



which is obtained from fits to the hadronic spectral function We empha- 
size, however, that there is abuse of notation in (^) since the field strength 
F^j, is defined in terms of the coset (gluon) fields which do not contribute to 
the above defined condensate O4 by construction. First results which employ 
the Wilson loop for extracting the gluon condensate were obtained in the 
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Figure 6: The gluon condensate O4 in units of the string tension a squared, 
as function of the coohng scale n/y/a (see text). 



early eighties from lattice gauge theories . The most recent value 

for the non-natural case of the pure SU (2) gauge theory is given by |^ 



O4 ^ 0.15 GeV^ (30) 

and was calculated from field strength correlation functions (for a recent 
review of the method see e.g. |3j 



In order to assign a physical meaning to the quantity (|28|), it is important to 



check whether this yields a finite and non-vanishing value in the continuum 
limit a — > 0. We have calculated /2 = 7r^G2a^/24 as function of f3 for cooled 
configurations, fixing the cooling scale = 0.5 cr. The result is shown in 
figure 1^. Our results nicely meets with the expectation from perturbative 
scaling (dashed line) for /3 > 2.2. This observation only tells us that we 
can safely extract a signal for the desired expectation value which remains 
valid in the continuum limit. Again, it turns out advantageous to have the 
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self-restriction of the cooling in a scale independent way. 

In order to show the stripping off of the gluonic contributions from the OPE 
parameter O4 (^), we can follow the limit k — 0. The result for O^Ik) ( ]28|) 
as function of n given in physical units is presented in figure |^. Data from 
various (3 in the scaling region (3 > 2.2 fall on the same curve. For n > ^/a, 
a significant contribution of the gluon radiation to the "gluon condensate" 
is still present, and one asymptotically expects O^i^K.) oc k,"^. This behavior is 
confirmed by our lattice data (see figure |^). Roughly at the confining scale, 
set by K ~ cr^^"^, the contribution of the gluons to Oi^n) becomes comparable 
with the gluon condensate carried by the Z{2) vortex vacuum texture. For 
clearly displaying the contribution of the texture, we compare the lattice 
data with the model fits 

= ao + ai k"^ , (fit A) 

OJa^ = h^K\ (fit B) . 

In both cases, the terms proportional to parameterize the gluonic contri- 
bution while the term of fit A specifies the Z2 vortex content. The lattice 
data clearly favorites fit A. We find it convincing that the "gluon conden- 
sate" Oi{K) approaches a finite value in the limit k — which is roughly 
consistent with the known value (|30D (dash-dotted line in figure |^). 

5 The positive plaquette model (PPM) 

The numerical results presented above were obtained with the SU (2) Wilson 
action, which includes a definite prescription of the interaction between center 
and coset fields. The question arises how the residual action density which 
is carried by the Z{2) gauge fields after 50(3) cooling does depend on the 
choice of the lattice action used in the simulation. In order to get some 
information on this dependence, we adopt an extreme point of view in this 
section and repeat the analysis of the previous section using the positive 
plaquette model (PPM). This model is defined by the partition function 

^ppm = J'^U H e{tr U^^^,) expj^tr f/,,^,} . (31) 

X flV 

The Gibbs weight is that of the Wilson action up to the fact that link config- 
urations which would lead to negative plaquettes are rejected. One expects 
that the latter constraint strongly influences the asymptotic Z{2) gauge field 
configurations remaining after cooling. 
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PPM model, 12 




Figure 7: The static quark anti-quark force calculated in the PPM. 



A set of low energy quantities - the string tension, the glueball masses and 
the topological susceptibility - have been studied in great detail within the 
PPM in [^|. It was found that such observables in units of the string tension 
become independent of the lattice regulator a when the continuum limit is 
approached within a certain scaling window (3 G [1.3, 2.1]. Moreover, the 
values of the above observables quantitatively agree with the values obtained 
with Wilson action. However, the renormalization group scaling of the lat- 
tice spacing with (3 does not match with the expectations from continuum 
perturbation theory in the investigated scaling window of the PPM. In fact, 
a scaling 



^a{(3) = 0.36 - 0.3 (/3 - 1.3; 



/?G[1.3, 2.1] 



(32) 

In order to test 



is consistent with the numerical data presented in 
our algorithm for simulations of the PPM, we have re-calculated the quark- 
antiquark force in physical units using a(/3) (P^) as input. We could verify 
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Figure 8: The SU{2) action density as function of the coohng scale k in the 
PPM. 

that the data points obtained for several values of the lattice spacing fall on 
top of the same curve (see fig. ^ . 

The continuum limit, however, is problematic to assess. Notice that we 
cannot assume that (|3^ ) is valid upto a = 0, otherwise the equation suggests 
a strong coupling UV fix point at 

^7fix := 9{a ^ 0) ^ 1.26 , {(3 = A/g^) (33) 

and a renormalization group /^-function of 



where we have defined the UV cutoff by A = 7r/a. Such a fix point contradicts 



asymptotic freedom. The equations (|33|j3^) are based on extrapolation of 
numerical data to the continuum limit. Therefore, it cannot be excluded 
that the onset of perturbative physics is postponed to the region P > 2.1. 
Note, however, that at /3 = 2.1 the UV cutoff is of order A 11 GeV (if we 
take the reference scale ^/a = 440 Me V), where perturbative scaling should 
dominate. 

Cooling (as well as gauge fixing) is a non-local procedure on the link con- 
figurations. Therefore, it is not excluded that we obtain negative plaquettes 
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after 5*0(3) cooling although we started from a configuration with positive 
plaquettes onlyQ. Formally repeating the analysis of the previous section for 
the PPM, we present our result for the "gluon condensate" as function of the 
cooling scale in figure ^ (left panel), compared with the corresponding result 
obtained with Wilson action in a double-logarithmic plot. Again the conden- 
sate value for each cooling scale does not depend on the bare inverse coupling 
/? of the PPM as long as the scaling (^) is applied. We find in the case of 
the PPM that the asymptotic value (for k — 0) is non-vanishing, but signif- 
icantly smaller than in the case of the Wilson action. In fact, we find that 
our data for O^Ik) are well reproduced by the fit (see the semi-logarithmic 
plot in figure ||, right panel) 

Oi{K) = 0.45^2 exp|2.63-^| . (35) 

Extrapolating (^) to k = 0, we find that the gluon condensate of the PPM is 
non-vanishing, but one order of magnitude smaller than the gluon condensate 
if the Wilson action is used. 

The fact that we have obtained a renormalization group invariant and non- 
vanishing gluon condensate also in the PPM, in which the plaquettes are 
constrained to be positive, stirs the hope that our above defined gluon con- 
densate gets physical significance independent of the choice of the lattice 
action. We attribute the fact that we don't find the same value of the con- 
densate in the case of Wilson action and in the PPM to the deficiency of the 
latter to match with perturbative scaling at very short distance. This is be- 
cause the gluon condensate appears as the first correction to a perturbative 
calculation. 

6 Conclusions 

In this paper, we have separated the SU{2) gauge field degrees of freedom 
into thin Z{2) center vortices and 50(3) coset fields. Since the 5*0(3) coset 
fields are isomorphic to algebra valued fields, these degrees of freedom have 
been identified with the gluonic ones. 

A new self-restricted cooling algorithm which reduces the 50(3) action of 
the coset fields facilitates the gradual removal of the gluon fields from the 
lattice configurations while preserving the center degrees of freedom. This 

^The same 5*0(3) action is used for the cooling, the positive plaquette constraint is 
ignored during the coohng. 
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SO{3) cooling procedure is gauge covariant. Hence, the remaining SU{2) 
Wilson action density reveals the gauge invariant Z{2) vortex texture of the 
SU{2) vacuum. 

Extracting the string tension for several values of the gauge invariant cooling 
scale k/x/o", we have found that the string tension vanishes in the limit of 
unlimited cooling. This shows that the Z{2) vortex ensembles remaining 
after SO {3) cooling cannot be identified with the confining P-vortices found 
in the MCG. 

The operator product expansion (OPE) does not offer an unambiguous pre- 
scription for identifying the condensates. We have suggested here an ap- 
pealing picture: the mass dimension four condensate is given by the action 
density of the effective Z{2) gauge model configurations which remain after 
the 5*0 (3) cooling procedure. 

This proposal gets support from the following numerical observations: first, 

the SU{2) action density for a given cooling scale properly scales towards the 
continuum limit; second, this action density approaches a rcnormalization 
group invariant constant in the limit of infinite cooling, when the SU{2) 
field is reduced to its Z{2) (vortex) content. This quantity gets immediate 
importance as the gluon condensate figuring in the OPE approach. 
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